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We show that according to the present understanding of the energy evolution of the observables 
measured in deep-inelastic scattering, the photon-proton scattering amplitude has to exhibit geo- 
metric scaling at each impact parameter. We suggest a way to test it experimentally at HERA. A 
qualitative analysis based on published data is presented and discussed. 



I. INTRODUCTION 

"on : 

The striking phenomenon of geometric scaling discovered in the HERA data has triggered a wide interest. The 
inclusive deep-inelastic scattering (DIS) cross section cr(x, Q) measured as a function of the photon virtuality Q and 
the Bjorken variable x, was shown to depend on the single combined variable Q 2 ■ (x/xq)^ where Ao ~ 0.3 and 
xq ~ 3 • 10~ 4 . Empirically, this scaling turns out to be valid for x < and for all available values of Q 2 (in practice, 
for Q 2 = to 400 GeV 2 ) . Geometric scaling has also been found more recently in deep-inelastic scattering off nuclei y] . 

A similar scaling behaviour was discussed a long time ago in the framework of elastic proton- proton scattering j^Qj. 
At that time, an impact parameter analysis revealed that the rise of the cross section as a power of the center-of-mass 
energy was due to an effective growth of the size of the interaction region, which was almost black in its center. 
(N ■ Geometric scaling in 7* — p scattering was discovered in the context of the discussion of saturation effects, and is 
often considered as an evidence for the fact that the proton looks already quite black, or "saturated" , at HERA. 

The scaling variable turns out to be the dimensionless ratio of the photon virtuality Q and the so-called saturation 
^vq ■ scale, which has the following dependence upon the Bjorken variable: 
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Q s {x) ~ 1 GeV • ( — ) (1) 
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The latter is believed to characterize the mean transverse momentum of partons with longitudinal momentum frac- 
O .' tion x. The growth of this typical momentum scale with energy W — Qj \fx is understood within several different 
formulations of dense parton evolution [H 0, IS El El El El which technically lead to non- linear integro-differential 
iq" equations. Geometric scaling was seen in numerical simulations of these models ll^JTiL It was shown to be consistent 
with QCD evolution in the range in which it is verified in the data In Ref.|lfit 1171 ITsj. it was derived from per- 
turbative QCD for Q above the saturation scale as a feature of the Balitsky-Fadin-Kuraev-Lipatov (BFKL) equation 
.£h ! (l9ll20j with appropriate boundary conditions. In Ref.[2l|. the symmetry of the quantity Q/Q s • cr(Q/Q s ) under the 
interchange of Q and Q s already noticed in Ref.Q] was interpreted as a possible evidence for a proton looking like a 
collection of independent dipoles of typical size 1/Q S . The generalized vector dominance model combined with the 
color dipole model as formulated in [23 also leads to geometric scaling. More recently, geometric scaling for all Q 2 
was deduced from the critical behaviour of the correlation function of Wilson lines which appears at small a; in a 
near-light-cone Hamiltonian formalism |23| . This approach looks fundamentally different from the other ones, in the 
sense that it does not rely on the parton model. 

Most of the latter discussions were only concerned with inclusive quantities. The proton was assumed an homoge- 
neous disk, and its size large with respect to the size of the probe, so that boundary effects could be neglected. The 
reason for these assumptions is the technical complication of the saturation equations when full dependence upon 
transverse coordinates is taken into account. 

However, one knows that the saturation scale itself depends on the impact parameter. Indeed, it is related to the 
density of partons inside the proton, which obviously decreases smoothly when one moves from the center to the 
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periphery. The saturation scale in the ideal case described before was only a kind of effective one. Actually, it is 
theoretically not clear what this scale is in a real proton, because all momenta from Aqcd to Q are a priori expected 
to contribute to cross sections. 

The imp ortance of impact parameter studies for the phenomenology of saturation effects at HERA was explained 
in Ref . |2J| . A theoretical discussion of the impact parameter dependence of the saturation scale Q s (x, b) was provided 
in In the latter paper, some arguments were put forward in favor of a "local" geometric scaling, which would 
manifest itself in the fact that cross sections for exclusive processes only depend on the ratio of scales Q/Q s (x, b). Let 
us also stress that geometric scaling for exclusive processes has recently been predicted qualitatively in the context of 
the discussion of generalized parton distributions j2|| . 

The aim of this paper is to show that "local" geometric scaling is testable experimentally. In the first section, we 
explain why this new scaling is expected. We then propose a method to investigate it at HERA. Finally, we discuss 
the qualitative results obtained so far by applying our method to the data. 

II. SATURATION AND GEOMETRIC SCALING 

In this section, we show how local scaling appears within different theoretical approaches. We consider the idealized 
process of scattering of a dipole of size 1/Q off a dipole of size 1/A. A and Q are fixed momentum scales, which are 
taken in the perturbative region in the following discussion, and the ordering Q > A is assumed. The small dipole is 
called "projectile" and the large one "target". In the end, A will be formally extrapolated to Aqcd- We denote by 
y = log(l/x) the maximum rapidity available at c.o.m. energy W = Q/y/x. 

A. Global and local scaling above the saturation scale from linear evolution 

We explore here the perturbative picture of high energy scattering. We consider dipole-dipole scattering at fixed 
impact parameter b as given by the BFKL equation. We assume that b is larger than the sizes of the initial-state 
dipoles: b^>l/Q, 1/A. The Mellin representation of the solution of the BFKL equation then reads [27ll28l| 

A d (y,Q,b) = \-QTj-^ J |^(l-2 7 )d( 7 )(l66 2 QA) 27 e a ^) , (2) 

where the integration goes over a line in the complex plane parallel to the imaginary axis and intersecting the 
real axis between and 1. In our conventions, the S and T matrices are related by S = 1 — T, so that generally 
speaking the amplitudes are essentially real at high energy. This is true in particular for the BFKL solution 
The following standard notations have been used: a s = a s N c /n and x( 7 ) = 2-0(1)— ^(7)— •0(1 — 7). The function 
d( 7 )=a^/(16 7 2 (l — 7) 2 ) is the Mellin transform of the elementary dipole-dipole elastic amplitude with respect to the 
ratio of their sizes, with appropriate normalization. 

The forward amplitude 7d is obtained by integrating (J2J over the impact parameter b. As formula J3J) holds only in 
a large b approximation, we must set the lower bound of the integral to a number of the order of the size of the largest 
interacting dipole. We choose 1/(4A) and check the consistency a posteriori. The forward amplitude then reads 

T A {y,Q)=2n r dbbA d {y,Q,b) = % [ ^(7) ( ? V ^ x(7) ■ (3) 
A/(4A) Q 2 J 2iir \AJ 

One checks the normalizations by computing the total cross section, using the optical theorem which reads a = 27^ in 
our conventions. Setting y = 0, one recovers the well-known expression of the elementary (lowest order) dipole-dipole 
cross section |29j . 

Let us recall how one obtains geometric scaling from this equation. One is interested in the large y region, where 
a saddle point gives the dominant contribution to the integral in Expanding the function x(l) to second order 
around 7 = 1/2, the amplitude 10) reads 

r d (,,Q) = ^^exp{^-Ilog ( QVA 2 )--i-log 2 ( gVA 2 )} , (4) 

with the definitions o; = 4o; s log2 (the BFKL intercept) and a 2 = 14£(3)a s . At fixed y, the amplitude becomes large 
(of the order of the area of the target) when the argument of the exponential is positive, which occurs below the scale 
Q = Q S (y) satisfying the relation 



Q 2 (y) = A 2 e Xy , with X = a 2 U/l + Auj/a 2 - l) 



(5) 



Note that this formula can be compared with the form of the saturation scale previously quoted provided one 
makes the identifications A — > Ao and A — > 1 GeV • x^ ^ 2 . The equation Q 2 — Q 2 {y) defines the so-called critical 
line below which, per definition, the amplitude (0 is large and non-linear corrections prove to be necessary to ensure 
probability conservation. The interesting feature of Eq. J3J is that it can be rewritten as a function of the ratio of the 
variables Q and Q s , up to some approximation. Putting Q and JSJ together, one arrives at 
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The factor which still contains the rapidity alone is slowly varying, hence it can be approximated by its value for a fixed 
typical rapidity tjq. Moreover, if Q is sufficiently close to Q s , one can trade the exponential for 1. This approximation 
applies when | logQ / 1 Q s \^a 2 y^Jl + Auj/a 2 . Introducing the notation 7= (-\/l+4u;/a 2 )/2, the amplitude rewrites 



which exhibits geometric scaling. 

To make contact with HERA phenomenology, one has to set A = Aqcd and to average the projectile size r with 
weight given by the squared photon wave function. This is the well-known dipole model [sol l3ll | which provides an 
appealing and useful picture of DIS in the leading-logarithmic approximation of QCD. The virtual photon-proton 
cross section reads 



<r 7 * P (y,Q) = / <fr^.( r ,Q) 8 ^ 7 .( r ,Q).27a(0,l/r) . (8) 

is the photon wave function on a qq colour dipole. Besides the transverse size r and virtuality of the photon 
Q, this wave function also depends on the longitudinal momentum fraction z of the antiquark. The notation "(g)" 
stands for an integration over z. For negligible quark masses, ^L(r, Q) <g> ijjj»(r, Q) is a function of rQ only which 
vanishes rapidly at large rQ . Thus the scaling property of the dipole-proton amplitude JJJ is transmitted to the 
photon-proton cross section itself which then exhibits the same behaviour as the amplitude 7d in Eq. Q . 

The formula just established would give a reasonable description of the data in a limited kinematical range, provided 
one considers u> and a 2 as effective parameters. Indeed, the values given by the leading order BFKL equation are 
known to be incompatible with the data on one hand, and to receive large corrections from next-to-leading order 
terms on the other hand. The parameters have to be adjusted such that A ~ Ao ~ 0.3 and 7 ~ 1 0, |33|. Then 
formulae (JIJ, J5j show that the choice xq = 3 ■ 1CP 4 implies Aqcd = 300 MeV, which is reasonable. To obtain more 
reliable expressions, neglected effects such as diffusion into the infrared have to be taken into account. As shown in 
Ref. > the latter have sizable effects both on the parametric form of the saturation scale JSJ) and on the expression 
of the amplitude Q), but do not spoil the overall picture. 

The previous analysis can be taken over to the impact parameter dependent amplitude J2J) in a straightforward 
way. The steepest-descent method gives (see |3flj ) 



Mv, Q, b) = 7 2 * /2 log(16^QA) exp <^ uy - log(166 2 QA) - — log^(166 2 QA) \ . (9) 
{■na z y) 6 l z y a A y J 

This amplitude is again of order 1 for Q = Q Sl the saturation scale Q s being given by 

Q ' (y ' &) = 25eW exp(Ay) - (10) 

Within the approximations already made in the inclusive case, the amplitude can then be expressed as a function of 
Q and Q s only: 

MQ,QAyM) = ^(2^)-(^-)\ (ii) 

where the saturation scale now depends on the impact parameter b. Thus we see that the amplitude for each impact 
parameter is a function of the ratio of the inverse size of the projectile and of the local saturation scale. 

However, the kinematical range in which formula may apply is quite restricted. Let us set A = Aqcd- We have 
required large impact parameters (in practice, b > 0.2 fm) and Q s close to Q, which has to sit in the perturbative 
regime (Q ^> Aqcd)- From formula IjlUII . one sees that the latter conditions require a large rapidity evolution, during 
which the saturation regime might already be reached for smaller impact parameters (i.e. Q < Q s {b') for some b' < b). 
It is not possible to quantify this effect without knowing the dynamics for all impact parameters. 

We shall now refer to another, maybe safer approach, which will however lead to comparable conclusions on "local" 
geometric scaling. 



B. The impact parameter profile as an initial condition 



In Ref.[25j, the argument that the quantum evolution was quasi-local in impact parameter was developed on the 
basis of an analysis of the non-linear evolution equations. As a result, the dipole-proton scattering amplitude at fixed 
impact parameter can be written in a factorized form, because the latter feature means that the impact parameter 
dependence remains unchanged throughout the rapidity evolution. The non-perturbative profile of the proton is 
introduced as an initial condition. The validity of such a factorization has been confirmed very recently in the context 
of 7* — 7* reactions [3] • The amplitude then reads 

A d (y,Q 1 b) = S(b)-T d (y,Q) , (12) 

where T d (y, Q) is the forward amplitude Above the saturation scale, T d (y,Q) exhibits geometric scaling as seen 
on Eq.0, thus the amplitude ^4d also scales. The profile function 

S(b) = -ml ■ e~ 2m - b (13) 

was postulated (we chose the normalization such that J d 2 b S(b) = l). It involves the long-distance scale l/rn w , which 
is put by hand and interpreted as an initial condition. 

This leads once again to a scaling of the form (|llfl . with a local saturation scale given by 

Q 2 s (y, b) = e - 2m " b / A • Ql{y) = e^ 6 /* . A 2 QCD e^ . (14) 

We see that this result is qualitatively close to Ea. HlOfl . except for the ^-dependence which is exponential rather than 
powerlike. This stems from the fact that non-perturbative physics with short range pion fields were introduced here, 
whereas previously only long-range Coulomb-like fields were considered. 

At this point, it may be useful to evaluate the saturation scale at HERA from these formulae. For values of the 
rapidity y of the order of those measured at HERA in the small- a; region (x^xq ^3 • 10 -4 ), one sees from formula i|14|) 
that Q s (b — 0) ~ 1 GeV. Q s reaches the value of Aqcd at the impact parameter 6~0.6 fm. This means that for a 
probe of size 1/Q<1 GeV -1 , the proton looks grey everywhere and extends up to a radius of about 0.6 fm. 



III. A WAY TO CHECK SCALING AT HERA 



As geometric scaling at fixed impact parameter seems to be a definite prediction of saturation models, it is worth 
testing it against the data. A former analysis [24[ has already shown the relevance of impact parameter dependent 
analysis in the discussion of saturation. There, the S'-matrix element for dipole-proton scattering at fixed impact pa- 
rameter was extracted from diffractive electro-production of vector mesons at HERA. This quantity can be interpreted 
as a transparency coefficient, and thus quantifies the "blackness" of the proton as seen by the projectile. Such kind 
of analysis is always model-dependent, since to get the correct normalization of the S'-matrix element, one has to rely 
on an ad hoc model for the final-state vector meson. Although it was shown in Ref. [24| that requiring that the initial 
state be a longitudinal photon limits the model-dependence, an estimated uncertainty of 20% was still recognized. 
However, dependence of the cross section on a scaling variable can be tested in a model-independent way. 

The scattering amplitude at fixed impact parameter can be extracted from differential cross sections for high energy 
quasi-elastic processes by a Fourier transform with respect to the momentum transfer. By quasi-elastic we mean that 
the initial and final states have the same number of particles. In practice, the process we are thinking of will be 
diffractive electroproduction of vector mesons. 

Let us consider such a quasi-elastic process. Its amplitude A is related to the differential cross section through the 
formula 

T t - < < 15 > 

where t is the usual Mandelstam variable and all other dependences are implied. At high energies, it is well-known 
that the scattering amplitude A is essentially real. Formula (jT5|) shows that A(t) is proportional to the square root 
of the differential cross section da/dt. A further Fourier transform with respect to the two-dimensional momentum 
transfer then gives the scattering amplitude for a fixed impact parameter 

f +oc ^- Fda 

A(b) = V^J o d(-t) Jo (by/% ^ , (16) 



where Jo is a Bessel function. For formula 1(16(1 to be valid, it is important that the quasi-elastic final state be well 
identified as a single particle or resonance. A more complex final state like "f*p — » ir + ir~p (which is what is actually 
observed in the experiments) can proceed through several different channels (with intermediate states p, w, /()•••)■ Its 
study would require the careful treatment of phases and interference effects, see for example Refs.[35|, |3(|. We shall 
assume that we have a single- channel process, and we use formula (|16|l to extract the 6-dependent amplitude from 
the available data. 

We shall now be more specific, specializing to the case of electroproduction of vector mesons j*p — + Vp, where V 
will in practice be a p meson. The picture is the following. At high energy, the photon splits into a dipole of fixed 
size r which scatters off the proton before recombining into a meson. Due to its high relative velocity, the dipole does 
not change size during the scattering. Accordingly, the scattering amplitude A is the convolution of the photon wave 
function ip 7 * , the meson wave-function ipy and the dipole amplitude Ad in the same way as in Eq.ljHJ) 

A(y,Q,b)= [ d 2 ri/l,(r,Q)®^ v (r)-A d (y,l/r,b) , (17) 



the forward amplitude 27^ appearing in (JHJ) being replaced by the impact parameter dependent one Ad- Due to the 
fact that the product of the wave functions ipL* ® tpv is peaked around a typical dipole size depending only on the 
external transverse momentum scales, the scaling of the dipole amplitude At] (see Eas-l jlll ) is transmitted to the 
photon-proton amplitude A. As it is well-known from both empirical tests |37ll38l| and experimental [39| studies, the 
relevant distance scale at the photon-meson vertex is a combination of the photon virtuality Q and of the meson mass 
M, namely 1/^Q 2 + M 2 instead of 1/Q as in the case of inclusive scattering. We then obtain from Eas. <jl()|) . (|14|) 
and HI 1(1 that the amplitude A in Eo. 1(17(1 should be a function of the scaling variable r = (Q 2 +M 2 )(x/xq) x only, once 
the impact parameter b is fixed. 

To check this statement, we apply the transformation 116(1 to the HERA published data on electroproduction of 
p° vector mesons. In the presently published experimental analysis the total cross section as well as the 

logarithmic t-slope B are quoted. The data can then be represented by the parametrization 

^(Q, W) = B(Q, W) ■ a(Q, W) ■ e ~ B ^™ , (18) 
at 

up to \t\ of the order of |£ max | = 1 GeV 2 , which corresponds to the range in which data have been collected. The 
scattering amplitude at each impact parameter reads 

A(y, Q,b) = ^U d(-t) J (bV\t\)VB^ ■ & - B ^' 2 + J d{-t) J {bV¥\)y ' ° 




(^ + °° d{-t) J (b^[T\)VB^ ■ e- s l*l/ 2 ^ = 2 



Ua{Q,W) „- b 2 /2 B 
B(Q,W) 



(19) 



■e 



/2B(Q,W) 



We take the total cross section from 40, 41] . The kincmatical range of these measurements is shown on Fig^ The 
experimental points range from Q 2 — 0.47 to 27 GeV 2 and W — 23.4 to 150 GcV. Due to the lack of data for the 
slope B for all values of Q 2 and W, we shall assume 

B(Q, W) = 0.6 • r _ ■ 9fi + 1 V 4a' Iog(W/75) , (20) 



(Q 2 + M 2 )°- : 

all quantities appearing in the previous formula being expressed in powers of 1 GeV. The parametrization 1(20(1 at 
W = 75 GeV is taken from jij. A logarithmic energy dependence has been added according to the Donnachie- 
Landshoff parametrization for the Pomeron 01 . We set the Regge slope a' to the value 0.25 which is measured in 
hadron-hadron cross sections. The relevance of this choice and its implications will be discussed in the next section. 

The results of our analysis are displayed on FigElfor three different impact parameters, b = 0.3, 1.0 and 1.5 fm. 
Needless to say, our plot is only illustrative of the method and, being theoreticians, we cannot provide a quantitative 
analysis. The error bars shown take into account only the error on the total cross section (AA/A= (Aer/er)/2) and the 
contour gives a rough idea of the total uncertainty, evaluated as if they were uncorrelated and small by the formula 

AA l*,.l*Bf»\ (21) 



A 2 a 2 B V B 



We assumed an error on the slope B of 5% which seems realistic regarding the quality of recently measured data (see 
the J/^> analysis of 01)- We see from ((21(1 that the uncertainty on B affects mostly the large impact parameter 
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FIG. 1: Kinematics of the data on the total cross section for diffractive electroproduction of p° mesons at HERA 40, 41]. The 
points represent the location of the data in the kinematical plane (Q 2 , W). Curves of equal scaling ratio r = (Q 2 + M 2 )(x/xo) x 
are shown. 



region. For small 6, an additional error would come from the fact that the data does not extend to arbitrarily large 
\t\, hence the exponential form (|19f) is not constrained for small b. This is the reason why we performed our analysis 
only down to impact parameters as large as 0.3 fm (see [24| for a more complete discussion). 

The results on Fig[3 arc consistent with "local" geometric scaling within uncertainties, for each of the displayed 
impact parameters. Although the kinematical range of the data is not very large, scaling is a non-trivial feature since 
data points of different Q and W overlap, as seen on Fig^ We stress that from the technical point of view, the 
fact that the scaling holds for various values of b comes from a subtle interplay between the (Q, W) dependence of a 
and B in formula Finally, let us note that the data points denoted by a cross correspond to very low values of 

Q 2 (0.47 GeV 2 ), and are kinematically quite separated from the other data. As they are measured with a different 
method, we do not make a point of the fact that they do not seem to lie on the same curve as the larger Q 2 data. 
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FIG. 2: Amplitude for 3 different impact parameters as a function of the scaling variable. The dotted contours indicate a rough 
estimate of the total uncertainty due to experimental error on both the total cross section a and the slope B as explained in 
the text, see formula The data points are derived from HI |4l| and ZEUS [it^ analysis of diffractive production of p° 

mesons. The error bars shown take into account the uncertainty on the measurement of the total cross section only. 



IV. DISCUSSION AND CONCLUSION 



In this paper, we have demonstrated that "local" geometric scaling, predicted by the saturation models for deep- 
inelastic scattering, is compatible with the HERA data. However, a consistent experimental analysis with careful 
treatment of uncertainties as well as an extension of the kinematical range is needed to draw a strong conclusion. In 
particular, an extension of the range in W would enable us to compare measurements sitting at very different Q 2 but 
at the same scaling variable r. This would constitute a definite test of this local scaling. 

Nevertheless, let us note the surprizing fact that scaling seems to work well for large impact parameters, very far 
from the center of the proton, and even outside the "grey" zone (b = 1.5 fm), where the concept of saturation scale 
makes little sense. This feature depends of course a lot on the form of the measured B, and particularly on its energy 
dependence, which was introduced in an ad hoc way in Ea. (|20|l . Interestingly enough, if the energy dependence 
were less steep for large momentum transfers y/\t\, scaling would also be slightly better verified at small impact 



parameters. Such a feature is likely to happen since it seems that the harder the momentum scales, the weaker the 
energy-de pen dence of the t-slope. This rule of thumb is most seen in large-i proton-dissociative photoproduction of 
p mesons [45], although of course it is not a quasi-elastic process and our analysis would not apply. Moreover, recent 
preliminary experimental results indicate that a' is smaller than the Donnachie-Landshoff value, and decreases with 
the scale Q 2 + M 2 . If scaling were confirmed also for large impact parameters, it might be an indication of some more 
general symmetry of high energy dynamics. 

One would in principle be able to plot all data for different 6 as a function of a unique scaling variable Q 2 /Q 2 (y, b), 
see formula (|11H . However, we checked that the parametrizations for the saturation scale proposed so far in Eqs.JSJ, 
(|10|l do not lead to a superposition of all points on a single curve. This means that the available approaches to the 
computation of the impact parameter dependence of the saturation scale are not good enough yet. 

The theoretical status of the impact parameter dependence of the saturation scale is still unclear, mainly due 
to technical difficulties in solving the full saturation equations. Better insights into these issues would help a lot 
the understanding of high energy scattering within QCD. Important issues like unitarity problems require a good 
understanding of the evolution in impact parameter [la 1.3 3| (see Ref. .46] for a model which predicts such an evolution). 
On the experimental side, diffractive vector meson production in DIS is being investigated intensively at HERA |47j . 
and more data will be available soon. We believe that it would be worth to plot the forthcoming data in the way we 
propose in this paper. 
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